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LONG  TERM  GOALS: 

Development  of  accurate  and  fast  advanced  statistical  and  dynamical  nonlinear  models  of  ocean 
surface  waves,  based  on  first  physical  principles,  which  will  improve  and  accelerate  both  long  term 
ocean  surface  waves  forecasts  and  prediction  of  strongly  coherent  events,  such  as  freak  waves,  tsunami 
and  wave -breakings. 

OBJECTIVES: 

Study  of  non-stationary  limited  fetch  waves  growth  in  presence  of  wind;  interpretation  of  altimeter 
data  through  study  of  self-similar  solutions  of  Hasselmann  equation;  studying  of  integrability  of  ID 
Zakharov  equation  for  surface  waves;  study  of  the  possibility  of  generalization  of  compact  ID  water 
waves  equation  for  2D  situation;  study  of  the  implications  of  modulational  instability  on  solitons, 
rogue  waves  and  air-surface  interaction;  study  of  ID  approximate  model  of  deep  water  surface  waves 

APPROACH 

Numerical  methods  for  solution  of  integro-differential  equations;  analytical  self-similar  solutions  for 
integro-differential  equations;  Hamiltonian  formalism;  comparison  of  analytical  and  numerical 
solutions  with  experimental  data;  analytical  and  numerical  solution  of  approximate  model  for  deep 
water  surface  waves 

WORK  COMPLETED: 

•  Non-stationary  fetch  limited  numerical  simulation:  classical  and  non-classical  regimes. 

•  Physical  model  of  sea  wave  period  from  altimeter  data 

•  Prove  of  non-integrability  of  1-D  Zakharov  equation  and  generalization  of  compact  equation  for 
almost  1-D  waves 

•  Study  of  modulational  instability  and  its  implications  for  solitons,  rogue  waves  and  air-surface 
interactions 

•  Approximate  model  for  ID  deep  water  surface  waves  based  on  MMT  equation 
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RESULTS 


1,  Non-stationary  fetch  limited  numerical  simulation.  Classical  and  non-classical  regimes. 

During  reported  year  we  eontinued  to  study  Hasselmann  equation  for  energy  spectral  density  s{r,  k,  t) 


ds  ^  5®,,  ds 


dt  dk  dr 


^ ill  ^ wind  ^ diss 


(1) 


for  the  simplified,  but  real  physical  situation,  which  includes  the  most  important  component  of  the 
general  statement  —  non-stationarity  in  time,  exact  nonlinearity  ,  wind  forcing  term  and  wave 

dissipation  term  due  to  wave -breaking. 

Such  simplified  formulation  is  the  limited  fetch  wind  growth  situation,  where  Hasselmann  equation  is 
reduced  to 


:^con0^  =  S„  +  S.,„+S„, 

2co  ox 

where  x  is  the  spatial  coordinate,  orthogonal  to  the  shore,  and  0  is  the  angle  between  individual 
wavenumber  k  and  the  axis  x. 

This  situation  is  presented  schematically  on  the  following  diagram; 


(2) 


Wind 


Fig.l  Schematic  description  of  the  energy  fluxes  along  the  fetch  in  real  and  Fourier  spaces. 
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Black  axis  represents  the  real  spaee  axis  x  ortogonal  to  the  shore  line,  red/green  eylinder  show  Fourier- 
spaee  of  the  system.  Red  part  of  the  eylinder  eorresponds  to  positive  adveetion  veloeity  in  front  of 
time  derivative  of  Hasselmann£'^.f2^  (red  arrow  is  direeted  away  from  the  shore),  while  green  part  of 
the  eylinder  corresponds  to  negative  adveetion  velocity  (green  arrow  direeted  toward  the  shore),  in 
eorrespondence  with  the  cos^  mEq.(2).  Sueh  schematie  pieture  suggests  that  limited  fetch  growth 
consists  of  essentially  three  different  proeesses  happening  in  real  and  Fourier  spaees:  (i)  energy 
adveetion  away  the  shore  in  real  spaee  (red  tube);  (ii)  energy  adveetion  toward  the  shore  in  real  space 
(green  tube);  (iii)  nonlinear  interaction  of  the  waves  between  red  and  green  parts  of  the  tube  in  Fourier 
spaee  for  any  given  point  of  the  feteh  and  time. 


For  the  solution  of  Eq.(2)  we  need  to  know  the  source  terms  in  its  right-hand  side.  The  proeedure  of 
exaet  caleulation  of  term  is  well-known  and  given  by  Webb-Resio-Traey  algorithm.  The 

parameterization  of  the  wind  input  term  is  given  by  (  see  [Rl]  )  : 


=0.159 
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m  = 


\  cos  6  for  -7iI2<9  <  nil 
0  otherwise 


(3) 


g  Pa 


=  1.3.10- 


u 
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and  is  based  on  the  faet  that  Eq.(2)  has  self-similar  solution 


with  parameters 


q  = 


P=hs=- 


(4) 

(5) 


along  with  the  faet  of  fitting  the  experimental  data  by  speeific  regression  line,  see  [R2],  [R3]. 


The  eontribution  of  the  dissipation  term  was  caleulated  similar  to  [R3],  where  white-eapping 


CO 


dissipation  term  was  introdueed  implicitly  through  /  ( /  =  —  )  energy  speetral  tail  stretehing  in 


In 


frequency  range  from  =1.1  to  =  2.0  . 


We  solved  Cauchy  problem  for  the  Eq.(2)  for  the  following  set  parameters 


•  Fetch  size  40  km 

•  Number  of  nodes  40 

•  72  logarithmic  frequencies 
•36  points  angle  resolution 

•  Wind  speed  U^Q  =5m/  sec 


which  translates  into: 
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•  Characteristic  angular  frequency  =  g  /  =  2rad  /  sec 

•  Characteristic  linear  frequency  /o  =  1 . 1/fe 

•  Dimensionless  fetch  y  =  e  =  1 .6  •  1 0"^ 

U 

'^10 

Initial  conditions  were  chosen  in  the  form  of  low-amplitude  noise  in  the  red  part  of  the  Fourier  space 
and  zero  in  the  green  part  of  the  phase  space  cylinder  (see  Fig.  1). 

Fig.2  represents  total  energy  of  the  fetch  as  a  function  of  time.  One  can  see  that  energy  growth 
evolution  can  be  split  into  two  parts  —  first  of  relatively  fast  growth  for  time  less  than  4  hr,  and  the 
second  part  of  relatively  slow  growth. 


Fig.2  Total  energy  of  the  fetch  as  a  function  of  time. 

Comparison  of  the  phase  space  distributions  of  the  spectrum  for  t=2hr  (  see  Fig.  5  a,b  ),  t=8hr  (see 
Fig.6  a,b)  and  t=24hr  (see  Fig.  7 a,b)  shows  qualitative  difference  in  their  shape,  too. 

For  earlier  times  t=2hr  the  shape  of  the  spectrum  has  the  shape  of  the  single  hump,  growing  with  the 
distance  from  the  shore  along  the  fetch.  For  later  times  t=8hr  and  t=24hr  the  shape  of  the  spectrum  is 
more  complex:  besides  classical  center  single-hump  of  the  energy  spectrum,  growing  away  from  the 
shore,  one  can  see  side  satellites  in  the  form  of  spikes,  corresponding  to  waves  propagating  along  the 
shore. 

Furthermore,  the  direction  of  those  waves  propagating  along  the  shore  is  getting  slanted  toward  the 
shore  with  the  distance  diminishing  toward  the  shore.  Closer  to  the  shore,  the  waves  are  propagating  at 
15  degrees  toward  the  shore  line.  This  observation  is  quite  remarkable:  it  means  that  long  enough 
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excitation  of  the  waves  in  time  by  the  wind  blowing  away  from  the  shoreline,  finally  excites  the  waves 
coming  against  the  wind  toward  the  shore. 

Fig,3  represents  distribution  of  energy  along  the  fetch  for  different  times.  One  can  see  that  energy 
behavior  for  time  less  than  4  hr  is  described  by  threshold-like  behavior  propagating  along  the  fetch,  in 
correspondence  with  [R4], 


Fig.  3  Energy  distribution  along  the  fetch.  Dotted  line  —  linear  function. 


At  time  around  4  hr,  the  energy  growth  is  approximated  by  linear  function  of  the  fetch  distance.  This 
dependence  was  predicted  in  [R5]  through  self-similarity  analysis:  it  was  shown  that  for  fetch-limited 


self-similar  solution  the  dimensionless  energy  E  = 


=  ax .  For  time  t  =  4  hr,  we  get  from 


numerical  simulation  results  a  =  2.4  •  10  ’ .  Table  1  represents  experimental  data  (see  [R5]  )  on 


measuring  the  constant  (I  for  different  experiments.  One  can  see  good  correspondence  of  our 
numerical  simulation  and  experimental  predictions: 
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Table  1.  Results  of  measurements  of  the  constant  tt  for  four  different  experimental  observations. 


Experiment 

a 

Nakata  Bay  (Mitsuyasu  et  al.,  1971) 

2.89  X  10-7 

JONSWAP  (Hasselmann  et  al.,  1973) 

1.6  X  10-7 

Lake  St.  Clair  (Donelan  et  al.,  1992) 

1.7  X  10-7 

Bothnian  Sea  (Kahma,  1981) 

2.6  X  10-7 

Fig.  4  shows  mean  frequency  distribution  along  the  fetch  for  different  moment  of  time  less  than  4  hr. 
One  can  see  that  the  portion  of  the  "threshold-like"  function  is  closely  described  by  self-similar 
soXvXion Eqs.(4)-(5)  : 


Fig.4  Logarithm  of  mean  frequency  distribution  as  a  function  of  the  logarithm  of fetch  for  four 
different  moments  of  time:  1  hr,  2  hr,  3  hr  and  4  hr.  Solid  line  —  self-similar  prediction  q=0.3 

from  Eqs.(4)-(5) 

All  the  above  facts  clearly  show  that  limited  fetch  growth  in  the  channels  of  final  width  can  be  split  in 
time  for  two  different  physical  regimes. 

The  first  regime  corresponds  to  self-similar  solution  developing  for  the  characteristic  times  defined  by 
ratio  of  the  channel  width  to  characteristic  advection  velocity  of  the  wave  field  and  could  be  called 
classical.  This  regime  corresponds  mainly  to  energy  advection  in  the  red  part  of  the  cylinder  on  Fig.l. 

The  second  regime  is  happening  later  in  time  when  the  energy  starts  to  reflect  from  the  opposite 
boundary  condition  and  corresponds  to  energy  advection  against  the  wind  in  the  green  cylinder  on 
Fig.l.  It  exhibits  itself  in  waves  generating  predominantly  parallel  to  the  shore  line,  tending  to  slant 
toward  the  shore,  as  approaching  to  the  beginning  of  the  fetch,  at  15  degrees. 
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To  wrap  the  results  of  researeh  up,  our  numerical  experiments  show  existence  of  2  regimes  of  wave 
turbulence  in  the  limited  fetch  conditions; 

•  classical  regime  as  a  combination  of  self-similar  regimes  for  moderate  duration  r  <  3  •  10"* 

•  non-classical  regime  exhibiting  in  waves  propagating  almost  parallel  to  the  coastal  line  (about 
15  degrees  near  the  shore)  analogous  to  Bose-condensation  for  r  >  3  •  10"^ 


Fig.  5a  Spectral  energy  distribution  as  a  function  of frequency  f  and  angle  6  at  the  fetch  distances  2 

km,  14  km,  26  km,  38  km  for  time  2hr. 
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Fig.6a  Spectral  energy  distribution  as  a  function  of frequency  f  and  angle  6  at  the  fetch  distances 

2  km,  14  km,  26  km,  38  km  for  time  8hr. 


tr\etq'^ 


Fig.  7a  Spectral  energy  distribution  as  a  function  of frequency  f  and  angle  6  at  the  fetch  distances 

2  km,  14  km,  26  km,  38  km  for  time  24.47  hr. 


Fig.  7b  Spectral  energy  distribution  as  a  function  of frequency  f  and  angle  6  in  polar  coordinates  at 
the  fetch  distances  2  km,  14  km,  26  km,  38  km  for  time  24.47  hr. 

Research  has  been  presented  in  [PI],  [P2]. 


2,  Physical  model  of  sea  wave  period  from  altimeter  data 

We  use  the  asymptotic  theory  of  wind  wave  growth  proposed  in  [R6]  that  predicts  the  Kolmogorov  - 
like  relationship  between  instant  total  wave  energy  E  and  total  net  wave  forcing  S  (wave  input  minus 
dissipation).  The  latter  can  be  associated  with  the  observed  rate  of  wave  energy  dE/dt.  Inversion  of  the 
so-called  weakly  turbulent  law  of  wind-wave  growth  (see  Eq.  (1.9)  in  [R6]  )  gives  a  simple  formula 
for  period  of  spectral  peak  Tp  as  a  function  of  significant  wave  height  //*  and  its  gradient.  The 
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relationship  can  be  used  for  processing  altimeter  data  assuming  the  wave  field  to  be  stationary  and 
spatially  inhomogeneous.  It  is  consistent  with  satellite  measurements  setup  and  is  written  as  follows 

(6) 

Self-similarity  parameter  in  Eq.(6)  is  set  constant  ass=0.67  as  found  in  recent  simulations  of  wind  wave 
growth  [R7],  Wave  height  derivative  along  the  wave  propagation  (group  velocity)  dHs/ds  enters  Eq.(6) 
in  very  low  power  1/1 0  that  reduces  dramatically  the  effect  of  the  unknown  wave  direction  relative  to 
the  satellite  track  where  the  derivative  can  be  estimated. 

In  contrast  to  all  the  models  of  wave  periods  mentioned  above,  Eq.(6)  does  not  contain  any  empirical 
coefficients  and,  thus,  does  not  require  calibration.  Additionally,  it  does  not  operate  with  backscatter 
coefficient  gq,  which  characterizes  sea  surface  roughness  in  gravity-capillary  range  and  cannot  be 
related  to  dominant  waves  straightforwardly. 

The  proposed  method  relies  upon  weakly  turbulent  mechanisms  when  instant  sea  state  is  determined 
by  external  flux  to/from  waves.  The  instant  state  and  its  apparent  perturbation  rate  recorded  by  spatial 
derivative  dHs/ds  allow  one  to  estimate  wave  period  in  certain  range  of  scales.  These  scales  should  be 
sulficiently  long  as  compared  to  the  relaxation  scales  of  weakly  nonlinear  surface  waves  and  relatively 
short  as  compared  to  typical  scales  of  variability  of  the  wave  field.  Typical  distance  between  two 
consecutive  counts  of  satellite  altimeter  5-7  kilometers  (1  second)  likely  satisfies  these  simple  criteria 
for  sea  waves  and  allows  Eq.(6)  to  capture  the  weakly  turbulent  mechanism  of  wind-wave  growth. 

The  new  algorithm  has  been  tested  for  the  data  collection  of  the  ESA  initiative  Globwave 
{http://www.globwave.org).  Different  empirical  models  have  been  compared  both  for  along-track 
records  and  in  terms  of  statistical  distributions  of  wave  characteristics  for  particular  ocean  regions. 


gT/(2. 


Eig.  8  Dependencies  of  non-dimensional  wave  height  on  non-dimensional  wave  period  (wave  age) 
calculated  by  different  methods  (see  legend)  for  800  km  track  of  JASON-2  in  the  Caspian  Sea, 

June  2,  2012. 
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Fig.8  presents  an  example  of  along-traek  evolution  of  non-dimensional  wave  height  and  period  (wave 
age)  ealeulated  within  different  approaehes  [R8],  [R9],  The  new  model  shows  good  agreement  with 
previously  proposed  empirieal  dependeneies.  Relatively  high  dispersion  of  the  new  method  (red 
points)  ean  be  partially  explained  by  intentionally  eoarse  data  proeessing  when  neither  smoothing  nor 
interpolation  of  data  has  been  used.  At  the  same  time,  strong  eollapsing  of  results  of  caleulations  with 
[R8],  [R9]  in  Fig.8  looks  suspieious  enough:  in  situ  experiments  do  show  rather  high  dispersion 
relatively  to  "referenee"  growth  eurves  like  one  of  [RIO].  As  our  study  shows,  the  "impoverishing" 
wave  dynamies  within  the  empirieal  parameterizations  of  wave  period  is  well  pronouneed  both  for 
along-traek  eases  like  one  presented  above  and  in  terms  of  statistieal  distributions  for  subsets  of  the 
Globwave  data  colleetion. 


This  study  presents  the  very  first  results  of  the  new  model  for  wave  periods  from  satellite  altimeter 
data.  The  model  is  free  of  any  empirieal  and  tuning  parameters,  therefore,  is  independent  on  features  of 
partieular  altimeter  deviees  (if  properly  ealibrated  in  wave  height  Hg)  and  regional  peeuliarities  (say, 
extreme  salinity  in  some  inland  basins). 

Agreement  with  previously  proposed  approaehes  points  to  the  validity  of  the  asymptotie  weakly 
turbulent  law  the  new  method  is  based  on.  While  having  very  few  in  situ  evidenees  of  validity  of  the 
weakly  turbulenee  theory  for  the  wind-driven  sea  for  the  last  50  years,  one  has  billions  satellite 
altimeter  reeords  that  are  eonsistent  with  conclusions  of  this  theory. 

Results  are  published  and  presented  in  [P3]-[P7]. 


3.  Prove  of  non-integrability  of  1-D  Zakharov  equation  and  generalization  of  compact 

equation  for  almost  1-D  waves 


We  studied  amplitudes  of  six-wave  interactions  for  compact  1-D  Zakharov  equation,  Hamiltonian  of 
which  is: 


H  =  ^b*d)j^bdx  -h  —  I  h'  I 


^(bb'*  -b*b')-k\b\ 


dx. 


where  b'  =  — ,  and  k  -  modulus  k  operator. 

dx 


It  was  found  that  six-wave  amplitude  (consisted  of  two  four-wave  amplitudes)  is  not  canceled  for  this 
equation  on  the  resonant  manifold.  Thus,  it  was  proven  that  1-D  Zakharov  equation  is  not  integrable. 

Also  we  have  presented  the  results  of  numerical  experiments  on  long  time  evolution  and  collisions  of 
breathers  (which  correspond  to  envelope  solitons  in  the  NLSE  approximation)  at  the  surface  of  deep 
ideal  fluid.  The  collisions  happen  to  be  nonelastic.  In  the  numerical  experiment  "nonelasticity"  can  be 
observed  only  after  many  acts  of  interactions  (collisions).  This  supports  "deep  water  nonintegrability" 
numerically.  One  can  see  small  radiation  after  100  collisions  of  two  breathers  (solitions)  on  Fig.  9: 
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Fig.9  Appearance  of  small  radiation  (red  line)  as  a  result 
of  100  collisions  of  two  breathers  (black  line) 

We  have  derived  generalization  of  eompaet  equation  for  almost  1-D  water  waves,  slightly  modulated 
in  transverse  direetion  with  the  Hamiltonian  as  following; 


H  =  ^  bdxdy  +  \\\K 


-ibb':-bX)-Kjb\ 


dxdy. 


Preliminary  numerical  experiments  on  freak-wave  formation  at  the  surface  of  3D  fluids  were 
performed  in  the  framework  of  this  model.  Typical  picture  of  the  surface  with  a  freak  wave  is  shown 
below  on  Fig.  10  ; 
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Fig.  10  Typical  picture  of  the  surface  with  a  freak  waves  as  a  result  of  numerical 
simulation  of  generalization  of  compact  equation  for  almost  1-D  water  waves 

Results  are  published  and  presented  in  [P8]-[P15], 

4,  Modulational  instability  and  its  implications  for  solitons,  rogue  waves  and  air-surface 
interactions 

We  studied  the  phenomenon  of  modulational  instability.  According  to  recent  studies,  the  modulational 
instability  seems  to  be  a  mechanism  that  is  responsible  for  a  number  of  phenomena  that  takes  place  in 
the  ocean.  In  particular,  experimental  and  numerical  works  have  shown  that  the  wave  breaking 
mechanism  and  the  formation  of  rogue  waves  may  be  caused  by  the  modulational  instability  process. 

In  [P16]  have  performed  direct  numerical  simulation  of  the  Navier-Stokes  equations  for  simulating  a 
two-phase  flow  (water  and  air)  to  study  the  dynamics  of  the  modulational  instability  of  free  surface 
waves  and  its  contribution  to  the  interaction  between  the  ocean  and  atmosphere. 

If  the  steepness  of  the  initial  wave  exceeds  a  threshold  value,  we  observed  wave -breaking  events  and 
the  formation  of  large-scale  dipole  structures  in  the  air.  Because  of  the  multiple  steepening  and 
breaking  of  the  waves  under  unstable  wave  packets,  a  train  of  dipoles  is  released  in  the  atmosphere; 
those  dipoles  propagate  at  a  height  comparable  with  the  wavelength.  The  amount  of  energy  dissipated 
by  the  breaker  in  water  and  air  is  considered,  and  contrary  to  expectations,  we  have  observed  that  the 
energy  dissipation  in  air  is  greater  than  that  in  water.  The  possible  implications  for  the  wave  modeling 
of  aerosols  and  gases  exchange  between  air  and  water  are  also  discussed  in  the  paper. 

In  [P17]  the  rogue  wave  solutions  (rational  multi-breathers)  of  the  nonlinear  Schrodinger  equation 
(NTS)  are  tested  in  numerical  simulations  of  weakly  nonlinear  and  fully  nonlinear  hydrodynamic 
equations.  A  higher  accuracy  of  wave  propagation  in  space  is  reached  using  the  modified  NLS 
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equation,  i.e.  Dysthe  equation.  This  numerieal  modeling  allowed  us  to  direetly  eompare  simulations 
with  reeent  results  of  laboratory  measurements  [Rll].  In  order  to  aehieve  even  higher  physieal 
aeeuraey,  we  employed  fully  nonlinear  simulations  of  potential  Euler  equations.  These  simulations 
provided  us  with  basie  eharaeteristies  of  long  time  evolution  of  rational  solutions  of  the  NLS  equation 
in  the  ease  of  near-breaking  eonditions.  The  analytie  NLS  solutions  are  found  to  deseribe  the  aetual 
wave  dynamies  of  steep  waves  reasonably  well. 

In  paper  [PI 8]  the  problem  of  existenee  of  stable  nonlinear  groups  of  gravity  waves  in  deep  water  is 
eonsidered  by  means  of  laboratory  and  numerieal  simulations  with  the  foeus  on  strongly  nonlinear 
waves.  Wave  groups  with  steepness  up  to  0.30  are  reprodueed  in  laboratory  experiments.  We  show  that 
the  groups  remain  stable  and  exhibit  neither  notieeable  radiation  nor  struetural  transformation  for  more 
than  60  wavelengths  or  about  15-30  group  lengths.  These  solitary  wave  patterns  differ  from  the 
eonventional  envelope  solitons,  as  only  a  few  individual  waves  are  eontained  in  the  group.  Very  good 
agreement  is  obtained  between  the  laboratory  results  and  numerieal  simulations  of  the  potential  Euler 
equations.  The  envelope  soliton  solution  of  the  nonlinear  Sehrodinger  equation  is  shown  to  be  a 
reasonable  first  approximation  for  speeifying  the  wave -maker  driving  signal.  The  short  intense 
envelope  solitons  possess  vertieal  asymmetry  similar  to  regular  Stokes  waves  with  the  same  frequeney 
and  erest  amplitude.  Nonlinearity  is  found  to  have  remarkably  stronger  effeet  on  the  speed  of  envelope 
solitons  in  eomparison  to  the  nonlinear  eorreetion  to  the  Stokes  wave  veloeity. 

In  paper  [PI 9]  we  show  experimentally  that  a  stable  wave  propagating  into  a  region  eharaeterized  by 
an  opposite  eurrent  may  beeome  modulationally  unstable.  Experiments  have  been  performed  in  two 
independent  wave  tank  faeilities;  both  of  them  are  equipped  with  a  wave -maker  and  a  pump  for 
generating  a  eurrent  propagating  in  the  opposite  direetion  with  respeet  to  the  waves.  The  experimental 
results  support  a  reeent  eonjeeture  based  on  a  eurrent-modified  nonlinear  Sehrodinger  equation  whieh 
establishes  that  rogue  waves  ean  be  triggered  by  a  non-homo geneous  eurrent  eharaeterized  by  a 
negative  horizontal  veloeity  gradient. 

The  paper  [P20]  is  a  review  paper  on  rogue  waves  in  different  physieal  aspeets  ranging  from 
oeeanography  to  nonlinear  opties. 

In  the  paper  [P21]  we  present  the  first  observation  of  dark  solitons  on  the  surfaee  of  water.  It  takes  the 
form  of  an  amplitude  drop  of  the  earrier  wave  whieh  does  not  ehange  shape  in  propagation.  The  shape 
and  width  of  the  soliton  depend  on  the  water  depth,  earrier  frequeney,  and  the  amplitude  of  the 
baekground  wave.  The  experimental  data  taken  in  a  water  tank  show  an  exeellent  agreement  with  the 
theory.  These  results  may  improve  our  understanding  of  the  nonlinear  dynamies  of  water  waves  at 
finite  depths. 

In  paper  [P22]  we  highlight  an  analogy  between  water  waves  and  eleetromagnetie  waves,  i.e.  the 
formation  of  a  supereontinuum:  starting  from  a  very  narrow  band  speetrum,  after  a  few  wave  lengths,  a 
very  broad  eontinuum  speetrum  (supereontinum)  is  observed  with  a  eonsequent  emission  of  solitons. 

5,  Approximate  model  for  ID  deep  water  surface  waves  based  on  MMT  equation 

We  studied  analytieally  and  numerieally  the  Majda-MeLaughlin-Tabak  (MMT)  model 

^  p,  “1^1  V^k  W +zt2  +*3 ^^1  ^^2 ^^3 
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adjusted  for  description  of  essentially  nonlinear  gravity  waves  on  the  surface  of  ideal  deep  water 
through  specific  choice  of  its  parameters  a  =  1/2,  P  =3  and  >.  =  1. 


The  water  surface  can  be  reconstructed  from  complex  function  y/{x,  t)  through  the  following  relation: 


nix,  0  =  ^  j  (y/^  +  y/\  )dk  (7) 

On  the  base  of  analytic  study  and  numerical  experiments,  we  observe  existence  of  robust  long-living 
objects  in  the  form  of  quasi-stationary  localized  structures,  which  are  the  quasisolitons. 

Quasisolitons  of  large  amplitude  turn  into  quasibreathers.  Their  amplitude  and  spectral  shape  oscillate 
in  time  (see  Fig.  11). 


Fig.ll  Dependence  of  the  soliton  maximum  max(|^(x,  t)f')  taken  over  integration  domain  [0,2;r] 
(  solid  line,  left  axis  )  and  the  second  moment  J(A:  -  dk  (  dotted  line,  right  axis  ),  on  time  t. 

[  k\y/i^^  dk 

The  average  wave-number  is  defined  as  k^  =  ^ - 

JKI  dk 

These  oscillations  at  the  time  moment  of  reaching  the  quasibreather  maximum  are  accompanied  by 
formation  of  weak  collapses,  exhibiting  their  selves  through  power-like  spectra,  which  can  be 
compared  with  “white  capping”  of  real  ocean  waves,  see  Fig.l2. 
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Fig.  12  Real  and  Fourier  space  distributions  of  wave  field.  Top  graph:  \y/{x,  as  a  function  of  x 

I  |2 

for  t=479.00.  Bottom  graph:  Fourier  spectrum  logjo|^^(x,  as  a  function  of  signed  logarithm  of 
wave  number  sign{k)\og^^}f\.  The  left  slope  of  the  spectrum  is  approximated  by  function  oc 
(dotted  line);  the  right  slope  is  approximated  by  function  oc  k  (dashed  line). 

Fig.13  presents  surfaee  elevations  for  the  same  time  as  Fig.12.  This  picture  looks  qualitatively 

similar  to  experimentally  observes  "Three  Sisters"  killer  wave  on  the  ocean  surface  [R12]  and  the 
resent  results  of  freakon  simulation  on  the  deep  water  surface  [R13].  These  slope  values  have  the 
meaning  of  the  original  Euler  equations  for  deep  waver  surface  gravity  waves. 
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Fig.l3  Surface  elevations  rjix,  t)  as  a  function  of  real  space  coordinate  x  for  time  t=479.00, 
corresponding  to  Fig.  12.  The  presented  data  are  dimensionless  and  scale  invariant  due  to  this 

property  of  the  corresponding  equations.  For  the  purposes  of  comparison  with  the  experiment,  one 

should  multiply  both  horizontal  and  vertical  axes  by  the  same  dimensional  factor,  corresponding  to 

particular  experimental  realization. 

The  research  has  been  published  in  [P23], 

REFERENCES 

[Rl]  V.  E.  Zakharov,  D.  Resio,  A.  Pushkarev,  New  wind  input  term  consistent  with  experimental, 
theoretical  and  numerical  considerations,  arXiv:  1212,1069  [physics,ao-ph] 

[R2]  Resio  D.,  Long  C.,  Equilibrium-range  constant  in  wind-generated  spectra.  Journal  of 
Geophysical  Research,  v.l09,  CO1018,  2004 

[R3]  Long,  C.E,  and  D.  Resio,  Wind  wave  spectral  observations  in  Currituck  Sound,  North  Carolina 
Journal  of  Geophysical  Research,  v.  1 12,  CO5001,  2007 

[R4]  Second  year  NOPP  ONR  report.  Nonlinearity  Role  in  Long-Term  Interaction  of  the  Ocean 
Gravity  Waves,  2012 

[R5]  Zakharov,  V.E.  Nonlinear  Processes  in  Geophysics,  12,  1011-1020,  2005 

[R6]  Badulin,  S.  I.,  A.  V.  Babanin,  D.  Resio,  and  V.  Zakharov,  Weakly  turbulent  laws  of  wind-wave 
growth,  J.  Fluid Mech.,  591,  339-378,  2007. 

[R7]  Gagnaire-Renou,  E.,  M.  Benoit,  and  S.  I.  Badulin,  On  weakly  turbulent  scaling  of  wind  sea  in 
simulations  of  fetch-limited  growth,  J.  Fluid  Mech.,  669,  178-213,  2011. 


20 


[R8]  Gommenginger,  C.  P.,  M.  A.  Srokosz,  P.  G.  Challenor,  and  P.  D.  Cotton,  Measuring  ocean 
wave  period  with  satellite  altimeters:  A  simple  empirical  model,  Geophys.  Res.Lett.,  30  (22), 
2150,  doi:10.1029/2003GL017743,  2003. 

[R9]  Quilfen,  Y.,  B.  Chapron,  and  M.  Serre,  Calibration/validation  of  an  altimeter  wave  period 

model  and  application  to  TOPEX/Poseidon  and  Jason-1  altimeters.  Marine  Geodesy,  27,  535- 
549,  2004. 

[RIO]  Toba,  Y.,  Local  balance  in  the  air-sea  boundary  processes.  Part  I.  On  the  growth  process  of 
wind  waves,  J.  Oceanogr.  Soc.  Japan,  28,  109-121,  1972. 

[Rll]  A.  Chabchoub,  N.  Akhmediev  &  N.  Hoffmann,  'Experimental  study  of  spatiotemporally 
localized  surface  gravity  water  waves'.  Physical  Review  E-Statistical,  Nonlinear  and  Soft 
Matter  Physics,  vol.  86,  no.  1,  pp.  016311  (2012) 

[R12]  C  .Kharif,  E.Pelinovsky,  A.Slunyaev,  Rogue  Waves  in  the  Ocean,  Springer,  ISBN  3540884181, 
2009 

[R13]  V.E.Zakharov,  A. I. Dyachenko,  About  shape  of  giant  breather,  Eur.J.Mech.B-Eluids,  29  (2010), 
127-131 

PUBLICATIONS 

[PI]  V.E.  Zakharov,  A.N.Pushkarev,  Nonlinear  generation  of  surface  waves  against  the  wind 
in  a  limited  fetch  growth  model .  20th  meeting  WISE.  April  21-25,  2013,  College  Park, 
Maryland,  USA  [oral  talk] 

[P2]  A.N.Pushkarev,  V.E.  Zakharov  Nonlinear  generation  of  surface  waves  against  the  wind 
in  a  limited  fetch  growth  model.  XXV  lUPAP  Conference  on  Computational  Physics, 

Moscow,  Russia,  August  20-24,  [oral  talk] 

[P3]  S.  I.  Badulin,  A  physical  model  of  sea  wave  period  from  altimeter  data.  Submitted  to  Journal 
of  Geophysical  Research  August  6,  2013. 

[P4]  S.I.  Badulin,  V.G.  Grigorieva,  V.E.  Zakharov,  On  relaxation  due  to  nonlinear  transfer  and 
physical  model  of  wave  periods  from  satellite  altimeters.  20th  meeting  WISE.  April  21-25, 
2013,  College  Park,  Maryland,  USA  [oral  talk] 

[P5]  Sergei  I.  Badulin,  Vladimir  E.  Zakharov,  Andrei  N.  Pushkarev,  Simulation  of  wind  wave 

growth  with  reference  source  functions.  Geophysical  Research  Abstracts,  Vol.  15,  EGU2013- 
5284-1,  2013,  EGU  General  Assembly  2013  [published] 

[P6]  S.  I.  Badulin  and  V.  G.  Grigorieva,  Weak  turbulence  theory  as  a  ground  for  measurements  of 
wind  wave  periods  by  satellite  methods.  International  Conference  ERONTIERS  OP 
NONLINEAR  PHYSICS  Nizhny  Novgorod,  Russia,  July  28  -  August  2,  2013,  [published] 

[P7]  S.I.  Badulin,  V.G.  Grigorieva,  A  Physical  Model  of  Sea  Wave  Period  from  Altimeter  Data, 
ESA  Living  Planet  Symposium,  9-13  September,  2013,  Edinburgh,  United  Kingdom 

[published] 

[P8]  A.I.  Dyachenko,  D.I.  Kachulin  and  V.E.  Zakharov,  "On  the  Nonintegrability  of  the  Pree 
Surface  Hydrodynamics",  JEEP  Letters,  V.  98(1),  (2013)  x>k(}).  43-47  [published] 


21 


[P9]  Dyachenko  A.I.,  Kachulin  D.I.  and  Zakharov  V.E.,  "Collisions  of  two  breathers  at  the  surfaee 
of  deep  water",  Natural  Hazards  and  Earth  System  Scienees,  Special  Issue;  Extreme  seas  and 
ship  operations,  (2013)  [publlished] 

[PIO]  Dyaehenko  A.I.  Zakharov  V.E.,  Shamin  R.V.  and  Badulin  S.I.,  "Ereak  waves  and  problems  of 
their  study".  In  the  book  "Mirovoi  Okean"  V.I.,  Ed.  E.I.  Eobkovskii,  p.  576-592,  2013 

[published] 

[Pll]  lUTAM  Symposium  2013:  Nonlinear  interfaeial  wave  phenomena  from  the  micro  to  the 
maero-seale,  Limassol,  Cyprus,  April  14-18,  2013,  A.I.  Dyaehenko,  V.E.  Zakharov,  D.I. 
Kaehulin,  "Compact  dynamical  equation  for  2D  water  waves"  [oral  talk] 

[P12]  International  Conference  "Landau  Days-2013",  Chemogolovka,  Russia,  June  24-26,  2013,  A.I. 
Dyachenko,  D.I.  Kachulin,  V.E.  Zakharov,  "Compact  equation  for  2D  water  waves"  [oral  talk] 

[P13]  V  International  Conferenee  "Erontiers  of  nonlinear  physics",  Nizhny  Novgorod,  Russia,  July 
28-August  2,  D.I.  Kachulin,  A.I.  Dyachenko,  V.E.  Zakharov,  "Analysis  of  integrability  of  the 
free  surfaces  hydrodynamies  equation  for  deep  water"  [oral  talk] 

[P14]  International  Conference  "ELUXES  AND  STRUCTURES  IN  ELUIDS",  St.  Petersburg,  June 
25-28,  2013,  Kachulin  D.I.,  Dyachenko  A.I.,  Zakharov  V.E.  "Nonintegrability  of 
Hydrodynamies  Equation  for  Deep  Water  Gravity  Waves"  [oral  talk] 

[P15]  XXV  lUPAP  Conferenee  on  Computational  Physics,  Moscow,  Russia,  August  20-24,  A.I. 
Dyachenko,  D.I.  Kaehulin,  V.E.  Zakharov,  "Ereak  waves  at  the  surface  of  deep  water"  [oral 
talk] 

[P16]  lafrati.  A;  Babanin,  A;  Onorato,  M  Physical  Review  Letters, Volume:  110  Issue:  18  Pages: 
184504  [published] 

[P17]  Slunyaev,  A.;  Pelinovsky,  E.;  Sergeeva,  A.;  Chabehoub,  A;  Hoffmann,  N;  Onorato,  M.; 
Akhmediev,  N  PHYSICAL  REVIEW  E  Volume:  88  Issue;  I  Artiele  012909  DOE 
I0.II03/PhysRevE.88.0I290  [published] 

[P18]  Slunyaev,  A.;  Clauss,  G.  E.;  Klein,  M.;  Onorato,  M.  PHYSICS  OP  ELUIDS  Volume:  25  Issue: 
6  Article  Number;  067105  DOE  10.1063/1.4811493  [published] 

[P19]  A.Toffoli,  T.  Waseda,  H.  Houtani,T.  Kinoshita,  K.  Collins,  D.  Proment,  and  M.  Onorato 
PHYSICAL  REVIEW  E  (Rapid  Comminication)  87,  051201(R)  (2013)  1539 

3755/20 1 3/87(5)/05 1 20 1  (4)  [published] 

[P20]  M. Onorato,  S.  Residori,  U.  Bortolozzo,  A.  Montina,  P.T.  Areechi  Physics  Reports  528  (2013) 
47/89  http ;//dx.doi.org/10.1016/j.physrep. 2013. 03. 001  [published] 

[P21]  A.  Chabehoub,  O.  Kimmoun,  H.  Branger,  N.  Hoffmann,  D.  Proment,  M.  Onorato  and  N. 
Akhmediev  PHYSICAL  REVIEW  LETTERS  Volume;  110  Issue;  12  Article  Number:  124101 
DOE  10.1 103/PhysRevLett.l  10.124101  [published] 

[P22]  Chabehoub,  A;  Hoffmann,  N;  Onorato,  M;  Genty,  G;  Dudley,  J  M;  Akhmediev,  N  Physieal 
Review  Letters  Volume:  111  Issue;  5  Pages:  054104  (Epub  2013  Aug  02)  [published] 

[P23]  A.Pushkarev,  V.E. Zakharov,  Quasibreathers  in  the  MMT  model,  Physica  D,  248,  55-61,  2013 

[published] 


22 


